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Abstract 

The transport properties of a random velocity field with Kolmogorov spectrum and time cor- 
relations defined along Lagrangian trajectories are analyzed. The analysis is carried on in the 
limit of short correlation times, as a perturbation theory in the ratio, scale by scale, of the eddy 
decay and turn-over time. Various quantities such as the Batchelor constant and the dimensionless 
constants entering the expression for particle relative and self-diffusion are given in terms of this 
ratio and of the Kolmogorov constant. Particular attention is paid to particles with finite inertia. 
The self-diffusion properties of a particle with Stokes time longer than the Kolmogorov time are 



determined, verifying on an analytical example the dimensional results of |Qlin.CD/0103018|. Ex- 
pressions for the fluid velocity Lagrangian correlations and correlation times along a solid particle 
trajectory, are provided in several parameter regimes, including the infinite Stokes time limit cor- 
responding to Eulerian correlations. The concentration fiuctuation spectrum and the non-ergodic 
properties of a suspension of heavy particles in a turbulent fiow, in the same regime, are analyzed. 
The concentration spectrum is predicted to obey, above the scale of eddies with lifetime equal to 
the Stokes time, a power law with universal —4/3 exponent, and to be otherwise independent of 
the nature of the turbulent fiow. A preference of the solid particle to lie in less energetic regions 
of the fiow is observed. 



PACS numbers: 47.27.Qb 47.55.Kf 02.50.Ey 



I. Introduction 

One of the differences between high Reynolds number turbulence and other examples of random 
fields with power law scaling, is the Lagrangian nature of time correlations . From the theoretical 
point of view, the need for a Lagrangian treatment of time correlations has been one of the main 
difficulties in the realization of statistical turbulent closures ||]. Because of this, many such theories 
assume from the start that the turbulence dynamics be equivalent to that of a random velocity field 
with identical energy spectrum but Eulerian time statistics, i.e. the fluctuations decay without 
being transported by the larger vortices ^, . Such an assumption does not work in the case of 
particle transport: both relative and self-diffusion are affected by the way in which time correlations 
are defined. 

Concerning self-diffusion, in Kolmogorov turbulence, fluctuations at a scale I within the inertial 
range, have characteristic velocity ~ li and decay time ^ l~i along fluid trajectories. Hence, in a 
time t the velocity of a fluid parcel will change by an amount of the order of that of a fluctuation 
with that lifetime, i.e. by . If the fluctuations were not advected by the flow, the fluid parcel 
would see the fluctuation only for the time ~ l^^ it takes to cross it. The variation of the fluid 
parcel velocity in a time t would be therefore ~ ti . 

Concerning relative diffusion, this process is determined by vortices with the size of the fluid 
parcel separation at the given time. If these vortices were fixed in space, their effect on relative 
diffusion would be proportional to the crossing time by the fiuid parcels, which is determined by 
the large scale properties of the flow. In other words, if time correlations were given in an Eulerian 
reference frame, the process of relative diffusion would not depend solely on the inter-particle 
distance and on the velocity difference, but also on the total velocity. 

Given the difficulty in defining a velocity field with Lagrangian statistics, a successful strategy 
for the treatment of transport has been to neglect time correlations altogether, i.e. to consider a 
velocity field u such that (uq(x, t)w/3(0, 0)) = C/Q/3(x)(5(t): the so called Kraichnan model [||. In 
this model, Eulerian and Lagrangian time statistics trivially coincide in what is the zero order of 
some perturbation theory in powers of the correlation time of the turbulence. It has been possible, 
in particular, to determine the anomalous scaling exponents of a passive scalar injected at large 
scales in the velocity field [Q, ^, ||, 0. The origin of this success is that, although the time structure 
of the velocity correlation is lost, that of the relative displacement, whose geometrical properties 
determine the passive scalar correlations, is preserved |l2[ |l3|. (For instance, particle pair 
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separation still obeys Richardson diffusion). 

The question, at this point, is how to introduce finite correlation times in a perturbative manner, 
but preserving the Lagrangian nature of correlations. There are practical reasons to do this. One 
motivation, of course, is to be able to determine the time correlations of the particle velocities. 
Lagrangian dispersion models |p^ , [T^ , ^ are based on the adoption of prescriptions on the form 
of these time correlations; to be able to determine them directly from the statistical properties of 
the velocity field would be therefore of some interest. 

It must be said that most of the prescriptions entering a Lagrangian dispersion model could 
be obtained, in practice, by dimensional reasoning or by experiments. In some cases, like in the 
presence of particles endowed with inertia, this turns out, however, to be a difficult task (l7[ |l8| . 
It is very difficult, for instance, to make assumptions on the preference of solid particles to lie in 



certain regions of the ffow instead of others |19, |20|, Solid particle transport by a turbulent 
ffow is an example of a situation in which careful treatment of the time dependent statistics of 
the velocity field is essential. It is precisely the interplay between the response time of the solid 
particle to the fluid, i.e. the Stokes time ts, and the characteristic times of the turbulent flow p2| , 
which determines the dynamics, and this is clearly lost when all the turbulent times are sent to 
zero. 

Recently, there has been strong theoretical interest on the problem of turbulence induced con- 
centration fluctuations in a heavy particle suspension. In ]2^ , the role of a finite correlation 
time of the turbulent field was recognized. In the case of particle with Stokes time in the 
turbulent viscous range was analyzed exploiting the fact that, in this case, the fluid velocity is 
spatially smooth on the scale of interest for the solid particle. In both |^ and |2^, however, the 
inertial range structure of the turbulent flow was disregarded altogether. The approach carried 
on here, allows instead to analyze the production of concentration fluctuations in any regimes of 
Stokes times, in particular in the inertial range, where qualitatively different behaviors for the 
concentration fluctuation build-up are observed. 

Purpose of this paper is to extend the Kraichnan model to short but finite correlation times, 
preserving, in a controlled perturbation theory, the Lagrangian structure of correlations, and pro- 
viding several applications to the transport of particles with and without inertia. The analysis will 
be confined to a situation of two-dimensional, stationary, homogeneous and isotropic turbulence. 

This paper is organized as follows. In section II, the equations determining the extension of the 
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Kraichnan model will be illustrated and their main properties discussed. Section III will be devoted 
to the dynamics of passive tracers; the self-diffusion and relative diffusion of fluid parcels, including 
the expression for the constants involved, will be determined; the effect of finite diffusivity will be 
discussed and the Batchelor constant for a passive scalar injected at large scale in the flow will be 
calculated. In section IV, the transport properties of a heavy particle with Stokes time longer than 
the Kolmogorov time will be studied, focusing on the relation between the correlation time for the 
fluid velocity sampled by the particle, and its Lagrangian and Eulerian counterparts. Section V 
will be devoted to calculation of the concentration fluctuations arising from compressibility of the 
heavy particle flow. In section VI, the bias introduced by inertia in the sampling of fluid velocity 
by solid particles (non-ergodic effects) will be analyzed. Section VII will be devoted to conclusions. 

II. Finite correlation time extension of the Kraichnan model 

A two-dimensional random velocity field with Eulerian correlation times scaling like the eddy 
turn-over time of a real turbulent flow can be obtained very simply, writing appropriate Langevin 
equations for the Fourier components of the vorticity field: 

dtq^it) + 7kqk{t) = hk^j,{t) (2.1) 

where qk{t) is the (space) Fourier transform of the vorticity: 

9(x,t)=Vi-u(x,i), (2.2) 

[for the generic vector v, we indicate v_l = (— W2,i'i)], and £,k{t) is the Fourier transform of a zero 
mean fully uncorrelated noise term of unitary amplitude: 

(C(x,i)^(0,0))=<5(x)<5(i) (2.3) 

The damping and forcing Kernels 7/5 and hk are chosen, for k <C r]~^, with r] the Kolmogorov 
length of the flow, as: 



while, for kr] > 1, some cut-off is imposed on the forcing amplitude Hk. In this way, the velocity 
spectrum Uk(i), defined by (uk(i)up(O)) = Uk(i)(27r)^ ^(k -|- p), wifl read, for krj <C 1: 

TTm-/L^<- 2k^exp(-7fe|f|) 

Uk(i) -4.C;,.,e3^^^^^-^ (2.6) 
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where CkoI and e play the role, respectively, of the Kolmogorov constant and the inertial range 
energy flux in a real turbulent field having this correlation spectrum. For fcg <C fc <C f?"^, we thus 
have the energy spectrum: — CkoI^^ k^s . Identifying 7^7 with the decay time and k^-Uf, ^ (0) 
with the turn-over time of an eddy at scale k~^, we see that p gives the ratio of the eddy turn-over 
and eddy decay time in the inertial range. The effect of sweep by the large scales, however, is not 
accounted for in this way. 

The most natural way to impose Lagrangian correlations in the random velocity field is to 
include an advection term in Eqn. (2.1), which will take the form in real space: 

(dt + u(x, t) ■ V)g(x, t)+J d2j/7(x - y)q{y, t) ^ j d^yh{^ - y)C(y, t) (2.7) 

This has the form of a vorticity equation in which the forcing and dissipation term, instead of 
being localized respectively at large and small scales, act over the whole of the inertial range, and 
this is reflected in their being nonlocal operators in real space. This is opposite to what happens 
in a real turbulent fleld, where energy balance is established between large scale forcing and small 
scale viscous dissipation, by means of the nonlinear cascade. A nonlinear cascade is still present 
because of the convection term, but it acts on the timescale of the eddy turn-over time, and, for 
large p, its effect is only a correction to that of the forcing and damping terms. Choosing p large 
has therefore the consequence that convection acts merely as a large scale sweep. 

Actually, Eqn. (2.7) looks a lot like the typical starting point of many turbulent closures 
^, ^ , in which 7^ gives the turbulent response function (eddy viscosity of small scales) and hk the 
nonlinear forcing by the cascade. For instance, p~^ coincides with the renormalizcd dimcnsionless 
coupling constant of the Renormalization Group (RNG) closure ||, |2^, and its smallness is there 
the basis for the establishment of a perturbation theory. Here, the philosophy is rather different: no 
parametrization of the turbulence cascade is sought, p is chosen arbitrarily large, and the similarity 
with real turbulence is expected to be only kinematic. (Also, the separation of a Kolmogorov 
constant out of the energy flux e is arbitrary). 

Things can be made a little bit more quantitative, introducing scale by scale the sweep time: 

n = ~ C^Je-iA|fc-i. (2.8) 

i.e. the time needed to a vortex of size k~^ to pass in front of a fixed probe. We see that sweep 
is important for all scales for which ^kTk < 1, i-c, from Eqns. (2.4) and (2.6), for k > kop^. The 
Kraichnan model is recovered when sweep can be neglected in all of the inertial range, i.e. for 
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p > (riko) 3 . This means basically that the zero correlation time limit is taken before the infinite 
Reynolds number limit ?//jo ^ 0. In this regime we have: 



V^it)^^^-^Cl^re^k-^Sit) (2.9) 



To understand what happens in the regime of dominant sweep, it is convenient to shift to La- 
grangian coordinates. Introduce then the coordinate z(i|x, to) of a fluid parcel which at time to is 
at X, and define the Lagrangian velocity: 

u^(x,t) = u(z(t|x,0),i) (2.10) 

and analogous expressions for q^{x,t) and the other fields. After introducing the increase of 
trajectory separation in a time t: Sz{t\x., y) = z(t|x, 0) — z(<|y, 0) — (x — y) Eqn. (2.7) becomes, in 
the new variables: 

dtq^i^,t)+ J d2y7(x-y + 5z(t|x,y))g^(y,t) = J dV(x - y + <5z(i|x, y))C(y, t) (2.11) 
which must be coupled with the equation for Sz; inverting Eqn. (2.2): 

dtSz{t\^,y) = d2r[G(x,r)-G(y,r)]q^(r,i) (2.12) 



with 



(x-r+^z(^Xll)k (2 131 



We see then that the natural expansion parameter of the theory is: 

feft|x,y) |u(x,0)-u(y,0)| ^1 

■ ■ ■ p 2.14 

|x-y| |x-y|7|x_y|-i 

i.e. the relative amount of particle separation increase in an eddy lifetime. The zero order of the 
theory, which is Gaussian and is described by Eqn. (2.6) after substituting u ^ u"", corresponds to 
neglecting trajectory separation in an eddy lifetime, while keeping the uniform large scale sweep, 
implicit in the Lagrangian field q^. 

Although the results which follow in the present paper are all obtained to the lowest order 
in the p expansion |26| , associated with neglecting all non-Gaussian effects in u, a diagrammatic 
expansion of Eqns. (2.12-13) in terms of the fields q^, Sz and their conjugate could be obtained 
by means of the Martin-Siggia-Rose formalism p7[ |. This expansion would only be valid locally 
around t = 0, since, at long times, trajectory separation becomes dominant. [To be consistent, this 
perturbation expansion should not receive contribution by correlations involving pairs of points 



in space-time such that 7|x-x'|-i 1^ ^ ^'1 > Pi but this is expected to be true from the exponential 
decay of the time correlations] . 

The interaction terms in the perturbation expansion are obtained Taylor expanding the kernels 
7, G and h (H working with the field action). The result for 7 is, for instance: 

7(x - y + Sz{t\x, y)) = 7(x - y) + ^ A^„7^"-'" (x - y)^^ (i|x, y)...fc.„ (t|x, y) (2.f 5) 

n=l 

with = 1 a coefficient which may scale when carrying on power counting. Similar coefficients 
Xg„ and Xh„ are introduced in the Taylor expansion for G and H . The theory is thus characterized 
by an infinite number of interactions involving vertices, which, to ©(p^"), have up to 2 + n legs. 
To check for divergences at large k in the perturbation expansion, we use power counting 



directly in Eqns. (2.11-13) |28 . Rescaling coordinates and times as: 

x^Ax and t ^ An, (2.16) 

Eqns. (2.9-11) remain invariant in form provided we rescale the various fields and interactions 
A = q^,Sz, A^„, Ag„, Ah„: as A Al^l^, with 

m - -I N = 1 

[A7J - [AhJ - [AgJ = 0. (2.17) 

This leads to expect logarithmic divergences at large k, meaning renormalizability of the field 
theory and the possibility of logarithmic correction to scaling, produced by renormalization of the 
parameters in Eqns. (2.11-13). 

It must be mentioned that marginal interactions and renormalizability are consequence of the 
dimensional relation implicit in Kolmogorov scaling: [q^] = — [t]. In general, had we set: 

7fe ~ r Hk^ (2.18) 

we would have obtained: 



2 ' ' 2 

I' 



[A7J = [AffJ = [AgJ = -nm + [9J) = -{s + 2-ir) (2.19) 



We thus see that super-renormalizability [A] < and non-renormalizability [A] > of the theory 
occur, respectively, for positive and negative [q^] + [t] |2^. This corresponds to the two regimes 
of eddy decay time becoming asymptotically longer (shorter) than the eddy turn-over time, and 
hence the nonlinearity becoming dominant (negligible) at large scales. 



Marginality of the interactions means that logarithmic divergences may arise both at large and 
small k. At small k, however, such divergences are not expected, due to the subtraction in the 
definition of 6z. The reason is sketched below (more details in a separate publication; it must be 
said, anyway, that this is not a surprise: Lagrangian closures |2] were introduced precisely to cure 
the infrared divergences arising in the original Eulerian theories). As it appears from Eqn. (2.17), 
small k divergence is due to internal lines in a loop diagrams involving the field Sz. The scaling of 
Eqn. (2.17) is associated with large, not with small k behaviors. In fact, the divergences occurring 
for large fc in a loop diagram will not change if we exchange (5z(t|x, y) z(i|x, 0)+z(t|y, 0) — (x+y); 
this because each small eddy contribute to the separation x — y an amount which is of the same 
order of the one to sweep. Now, the logarithmic divergence predicted at small fc in a loop diagram 
comes indeed, from equating the scaling of the sweep z(t|x, 0) + z(t|y,0) — (x + y) with that of 
the trajectory separation (5z(i|x,y) also at small fc, which is incorrect. For small fc, this scaling 
should be corrected by a factor fc per field 6z involved in the lines of the loop, and this is enough 
to eliminate divergence. 

III. Passive tracer transport 

1. Self-difFusion of a fluid parcel 

Lagrangian correlation functions in the form (u^(x, t)u^(x, 0)) = (u(z(t|x, 0), t)u(x, 0)) are the 
simplest objects one may try to calculate from the random velocity field introduced in section II. 
The starting point, to lowest order in and after sending the Kolmogorov scale 77 to zero, is 
the following modification of Eqn. (2.6): 

TT^m A^r .1 ^J-^^ exp(-7fc|t|) 

U,(t) = 4.C^„..3_^___^ (3.1) 

The Lagrangian correlation time is then readily calculated: 

rE' = (|u^n[ J di(u^(x,i).u^(x,0))]"' =2pcLe-U| (3.2) 

and we have the following relation between the turbulence level — {v?) — (|u^P) and the 
integral scales of the flow fco and t^: 

ul = SCxoieh^^ = dpcLi^TL (3.3) 

The correlation time tl is determined by the particular form of UJ^ we have chosen at small fc, 
which is non-universal. It is more interesting, and relevant from the point of view of Lagrangian 
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dispersion modeling [ p4| , p9[ , to calculate the Lagrangian time structure function: 

(K(x,<) -7.^(x,0)][w;^(x,t) - w;3(x,0)]) = i(|u^(x,t) -u^(x,0)|2)<5„^ (3.4) 

We discover immediately that, in order to have a self similar spectrum for the inertial range, the 
time correlations should have continuous time derivative at i = 0, a property not satisfied by Eqn. 
(3.1). 

This self-similarity violation can be illustrated in a simple way imagining the turbulence field 
in the neighborhood of the fluid parcel as a superposition of nested eddies with scale In, velocity 
Un and eddy turn-over time r„: 

/„ = /o2"", Un = uo2-^, r„ = To2"^ (3.5) 

If the time correlation decayed linearly for i ^ 0, we would have: 

(|u^(x,i) - u^(x,0)p) ^Y.ul-+Y,ul^ ul\og{To/t)-; (3.6) 

'^n ^ . 'HD 

Tn<t Tn>t 

Thus, identical scaling of and t„, and linear decay of correlations cause the largest space scale 
to contribute to the structure function at arbitrary short time separation t, in the same way as a 
vortex with eddy turn-over time t„ ^ t, whence the logarithmic correction involving tq. 

In order to have a quadratic behavior of the time correlation at i = 0, it is necessary that the 
noise ^ in Eqn. (2.7) be correlated in time, and the correlation must again be given along the 
trajectories. The appropriate modification to Eqn. (2.7) is therefore: 



{dt + u(x, t) ■ V)q(x, t)+J d2t/7(x - y)(7(y, t) = r(x, t) 

{dt + u(x, t) ■ V)r(x, t)+J d2y7(x - y)r(y, t) = J d^yh{x ~ y)e(y, t) 



(3.7) 



where, for k rj ^: 

7fc=pcLe^fc^ and H, = \hk\^ = 87Tpp{p + p)Cl^iel {k^ + k^f^ (3.8) 

It is easy to show that also the field theory associated with Eqn. (3.7) is characterized by marginal 
interactions: [A-y^] — [A^„] — [Xh„] — [Ag„] = and the considerations in section II extend to the 
present case. 

The zero order of the theory leads to the following correlation function: 

Ui^W-'^"^ 4^^^°^^-^/^-"'^'-^^-"'*' (3.9) 

and the time correlation has a quadratic maximum at i = 0. Calculation of the Lagrangian 
correlation time leads to the same result of Eqn. (3.2), with the substitution p while 



smoothness of the time correlation eliminates the logarithmic correction to the scaling of the 
Lagrangian time structure function. This structure function obeys in fact, after sending kg 0, 
the expected normal diffusion behavior: 

(|u^(x,<)-u^(x,0)|2) = 2Cof-|t| with Co = c|„,-:^logp/p; (3.10) 

p — p 

the constant Co is 0{p) and, as expected from the discussion leading to Eqn. (3.6), diverges 
logarithmically for p/ p oo. 

2. Relative diffusion 

Analyzing the transport of a cluster of particles requires consideration of time intervals, during 
which the space separations involved cannot be approximated as constant. Over these timescales, 
the short correlation time limit leads to a perturbation scheme, which treats the velocity field to 
zero order as a white noise. 

We focus on the case of a pair of particles. We have to study an equation in the form: 

dt[z^it\ro, 0) - z„(t|0, 0)] = ul{ro,t) - u^(0, t) = C/,/3(z(t|ro, 0) - z{t\0, Ompit) (3.11) 

with (^a(^)C/3(0)) = ^afjSit) and Uap to be determined. Due to the multiplicative noise nature 
of this equation, attention must be paid to the possible presence of drift terms arising from the 



Stratonovich prescription implicit in its definition 30 . It is easy to show that this drift is identically 



zero, either by direct calculation of the increment (5z(t|x, 0) for t in the inertial range, or noticing 
that: 

(5z(t|ro, 0)) - r dT[((u^(ro, r)) - (u^(0, r))] = 0; (3.12) 

JQ 

this, because of homogeneity of turbulence. For this reason, the separation process is described 
simply by: 

dt{M\ro, 0) - z„(t|0, 0)][zp{t\ro,0) - zp{t\0, 0)]) = D^0{z{t\ro,O) - z{t\0, 0)) (3.13) 

with 



D^pir)^ y dt(K(r,i)-<(0,i))][u;^(r,0) -7.^(0,0))]) (3.14) 

This tensor is easily calculated from Dii{r) for r = (r, 0), exploiting incomprcssibility. Using 
d0 sin^ 6* sin^ (x cos 6*) — ^[1 ~ Jo(2x) — J2(2x)], we find, in the limit fco ~^ 0: 

1 _i 

Dn(r) = ' (3.15) 
P 
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where 

dx x-i[l~~ 3o{x) - h{x)] ^ 0.265, (3.16) 





with J„ the Bessel function of the first kind, is evaluated in terms of Gamma functions [3l| ] using 
the formula dx x'^J^lx) — 2'' y^^t-|y|^t^ • From incompressibility we find therefore: 

1 _i 

[^+3i'5"/3-^jJ (3.17) 

We want to study the asymptotics of the separation process of two particles in the inertial range. 
The procedure is standard (see e.g 10); we introduce the distribution P for the separation r at 
time t, which will obey the diffusion equation (the summation over repeated indices convention is 
adopted throughout the paper): 

dtP = ^da^dpD^l^P (3.18) 
and look for an isotropic similarity solution in the form 

P(r, t) = t-^f{t-^r) = t-^f{R) (3.19) 

Equation (3.18) takes then the form 

1 _i 

MRc.f)+ \ d^mRM = (3.20) 
2 2p 

This equation has an unphysical solution, which is divergent in i? = 0, and a finite one: 

2 

/(i^)^exp(- ] (3.21) 

whose moments are: 

1 1 

(i?") = i?i+"di?/(i?) = I ( -^gj"' ) r(3 + ^) (3.22) 

with r the standard gamma function. 

From here, the expression for the particle space separation is obtained in a straightforward 
manner; for 72^-1^ ^ 1, indicating r{t) — z(t|ro,0) — z(t|0,0): 

10240a3^c|„, 

{r'{t)) = ce-^^ c = (3.23) 

i.e. the space separation obeys Richardson diffusion. For the relative velocity, we have, from Eqn. 
(2.6): 

(K(ro, 0) - <(0, 0)]2) = 2a| C^oic^ (r^ (t)) (3.24) 
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where 

poo 

dx x^i [1 - Jo(a;) - J2(x)] ~ 2.149 (3.25) 





and, using Eqn. (3.22), for jx-'^t ^ Ij we find the normal difi'usion behavior: 

3 

(K(ro,0)-<(0,0)]2)=ce-i, c= %^ (3.26) 

Passing to the smoothed out in time version of the velocity field provided by Eqn. (3.7), is 
accomplished, as in the case of tl, by exchanging p J+J- -'^^ lEH' ^'^^^ ^ sub-exponential 
behavior for the function f{R) and Richardson diffusion were observed in a direct DNS (direct 
numerical simulation) of two-dimensional turbulence in the inverse cascade regime. Based on the 
results of that paper, extrapolating applicability of our leading order expressions in p would give 
then (taking also p — > oo): p ~ 2. 

3. The role of difFusivity and the Batchelor constant 

The dynamics of passive tracers, contrary to that of fluid elements, feels the effect of molecular 
diffusivity. Due to finiteness of the turbulent correlation times, this effect does not consist purely 
of an additive noise contribution to the tracer velocity. Indicating by a the molecular diffusivity, 
the passive tracer velocity will have the form: 

v(x,t) + (2r7)^|(x,<) (3.27) 

with (^Q,(x, t)^^(0, 0)) — 6a(3S{'K)5{t) and v obeying an equation in the form: 

(9t+v(x,t)-V)q„(x,t)-|- J d^y-f{x-y)q,{y,t)~ J d^yhi^ - y)Ciy,t) 

= -(2a)5(^(x,t)- Vq„(x,i))4 ~CTV2(7„(x,t) (3.28) 

where qv — • v, (.)^ is an average limited to the noise ^ and use has been made, in converting 
the advection by molecular noise into a diffusion term, of Ito's lemma |Q. We see (it is assumed 
that the limit 77 ^ is already taken) that there is a renormalization of the damping kernel 7: 

7fe ^ 7fc + f^fc^ (3.29) 

which leads to a cut-off for the velocity at the inverse diffusive scale 

V^' = ipcLf^-eia-i (3.30) 
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We have then: 

(vk(t)v_k(0)) = 2^Piz^^(uk(i)u_k(0)), (3.31) 

and for small space separations r/rja- 0, we have a quadratic behavior for the velocity structure 
function: 



Jo 



^\cKoieh.~'r^\logr/v.\ (3.32) 
The transport of a passive scalar ^(x, t) will be described by the equation 

{dt + v(x, t) ■ V)0(x, t) = aV^0{x, t) + /(x, t) (3.33) 

with /(x, t) a source term. An interesting quantity to calculate is the fluctuation spectrum for 6 
in the case / is random in time and concentrated at large scale: 

(/(x + r,t)/(x,0))=F(r)5(0, = (f'' I'" (3-34) 

10, kor > 

We can thus consider {9) = 0. The equation for the steady state passive scalar correlation 0(r) = 
{e{x + r, t)e{x, t)) will then be, for kor <C 1: 

(6'(x, f)[v(x + r,t)- v(x, t)] ■ V6'(x + r, t)) = 2aV^e{r) + 4ee (3.35) 

For r — > 0, the left hand side of this equation is zero; we thus obtain ee = f (|V^|^), i.e. ee is 
the dissipation of passive scalar fluctuations. Following the same approach of the previous section, 
the velocity difference v(x + r,t) — v(x, t) is approximated by a white noise. From Ito's Lemma, 

its contribution in Eqn. (3.35) will be an eddy diffusivity £)^^(r), whose expression will coincide, 
for r ^ r^cr^, with the one for Dap provided by Eqns. (3.16-17). Drift terms coming from the 
Stratonovich prescriptions are ruled out with the same arguments used in the previous section. 
The resulting diffusion equation will then read: 

dadfsi^Dlpir) + 2adai})eir) + 4ee = (3.36) 

For r ^ r]cr D^- is essentially a correction to the molecular diffusivity, and will read, from Eqn. 
(3.31): 

Dl{v)= j dt(K(r,f)-t;„(0,t))][t;;3(r,0)-t;^(0,0))]) 



Tar 13 , 13 



+ y(^«/5-^)J^(^/'?-)^ (3-37) 
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For ?7ct -C r, D^j is approximated by Eqn. (3.17), the molecular diffusivity a can be neglected and 
Eqn. (3.35) takes the form: 



e" + T^e' = ^ (3.38) 

Solution of this equation gives automatically the passive scalar structure function ([6'(x + r,t) ~ 
6'(x, t)]^) = 2[6(0) — 0(7')] in the inertial range for 8: rja <^ r <^ k^^ . This structure function 
scales like ri and can be written in the form: 

([^?(x + r,i)-0(x,i)]2) = -?f^ (3.39) 

with the parameter B — ^ the so called Batchelor constant of the flow. As with relative diffusion, 
the case of a velocity field with smooth time correlation described by Eqn. (3.7) is recovered 
substituting p with . 

IV. Solid tracers: 1-particle statistics 

We consider the simplest case of a linear drag. In the presence of gravity (or of a constant 
external force) and of the turbulent velocity field u(x, t), the solid particle coordinate z^(t|x, 0) 
will obey the equation of motion: 

z^(t|x,i) =v^(x,i)+UG, z^(0|x,0)=x (4.1) 

where uq is the gravitational drift, that we suppose constant and uniform and is the fluctuation 
in the Lagrangian solid particle velocity, which obeys the linear relaxation equation: 

v-(x,i) =r5i(u(z-(t|x,0),i) - v-(x,<)) =r5i(u^(x,t) -v-(x,i)), (4.2) 

with Ts the Stokes time. (For a spherical particle of radius a and density pp, in a fluid of density pa 
and kinematic viscosity v, we would have: ^7 11 ^ Pp/ Po\] we are disregarding any effect from finite 
particle Reynolds number [^). From now on we shall identify Lagrangian quantities calculated 
on solid particle trajectories by the superscript P. 

In general the non-coincidence of fluid and solid particle trajectories makes the analysis of Eqns. 
(4.1-2) a very difficult task. The short correlation time limit p 00, however, allows to proceed 
perturbatively in the fluctuating part of the trajectory separation u^t + z(t|x, 0) — z^{t\x, 0). The 
physical motivation for this is that, from Eqn. (4.2), u^i + z(i|x, 0) — z^(t|x, 0) fluctuates on 
timescale ts with velocity scale fixed by those eddies which have decay time ts- Hence, for p 
large, the fluctuating part of trajectory separation remains small on the scale of these eddies. 
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Furthermore, when either uct > 6z{t\x + u^t, x), or 7|„Qt|-ii > 1) in other words, when either 

3 3 

^g' < 1 or ^g' > p""^ (provided p > 1, one of the two conditions is always satisfied), it is 
possible to approximate z(t|x, 0) + uct — z{t\x + uat, 0). 
To lowest order we have therefore: 

u(z^(i|x,0),i) =u(uGt + z(i|x,0),i) =u^(x + UGt,t) (4.3) 

We obtain immediately the fluctuation amplitude of the velocity difference between solid and fluid 
particle at a given position. From Eqns. (4.2-3) we can write: 

V = / 7^ / — "k(^) exp( + ik • x) (4.4) 

and from here we obtain, using Eqn. (3.2): 



where 



{{Va - Ua){Vi3 - U/})) = 6al3ul — „ i > a/3ul log{TL / Ts) (4.5) 

Jl X(l + ^^Xs) ^s'^-'-L 

^ TT. ' 



Us 



m'ur (4.6) 



for TS < tl, is the velocity scale of eddies with lifetime ts and ut is the turbulent velocity defined in 
Eqn. (3.3). In order to proceed to next order, it is necessary to calculate the trajectory separation: 

z^(t|x, 0) - z(t|x, 0) = uct + (1 - e"*/^^) / dr e^/^^u^(x, t) 

J — OO 

- /*dTCxp(-^^)u^(x,r) (4.7) 

Jo TS 

We notice from this equation that the inertia produced part of trajectory separation does not grow 
indefinitely. In other words, if ug = and to lowest order in p~^, there will be localization of 
solid particle trajectories around the fiuid parcel trajectories they cross at any given time; from 
Eqn. (4.7): (|z'°(t|x, -oo) - z(t|x, -oo)|^) ~ {uttsY ~ Cxo/e^ fco'^rj. We thus introduce the 
localization length Si: 

Sl=Cl,ieh-Ks (4.8) 

What happens is that the velocity difference v'' — obeys a relaxation equation with a forcing 
which is a time derivative; from Eqn. (4.2): ^(v'' — u^) + t^^(v^ — u^) = — u^. The frequency 
spectrum of v'° — does not have therefore the small frequency singularity necessary for long time 
divergence. The localization length Si will appear to play a fundamental role in the production 
both of concentration fluctuations and of corrections to the velocity correlation time. (Of course, 
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to higher order in p^^, the relative separation of fluid parcels sets in and localization is destroyed; 
Si becomes then, that part of trajectory separation which remains after the Richardson diffusion 
contribution is subtracted out). 

In the absence of gravity, beside the integral scale dependent localization length Si , three more 
scales, which, if <C r^, are purely inertial, can be obtained combining r^, the crossing time of 
an eddy by a solid particle, the eddy lifetime and the eddy turn-over time. We have the size S of 
an eddy whose lifetime equals Tg: 75-1T5 ^ 1; the size Sc of an eddy that is crossed by a solid 
particle in a time ts'. us ^ Sc/ts; the size Si of an eddy whose lifetime equals the crossing time 
by a solid particle: 5*^7^-1 ~ U5. Summarizing: 

S^p'^cioi-^^rl, S^^p^Cl^fe^rl and 5. = p-icLe^ (4.9) 

From Eqn. (4.9), we identify the following sequence of ranges: 

• A large separation range r > S, in which the fluid velocity u'' varies slowly on the scale of 
the relaxation time ts- 

• A first intermediate range 5 < r < S'c in which the fluid velocity u'' is a fast variable, but 
still, Ts is short compared with the crossing time of an eddy of size r; hence, Eqn. (4.2) has 
the form of a Langevin equation with a noise t^^u'' of constant amplitude on the scale of 
this crossing time. The cross-over scale S will play an important role in the determination 
of the degree of non-ergodicity of the solid particle flow (see section VI) . 

• A second intermediate range Sc < r < Si, in which the crossing time is shorter than both 
the Stokes time and the eddy turn-over time, but is longer than the lifetime of an eddy of 
that size; hence, the solid particle moves ballistically with respect to the fluid. 

• A small separation range r < Si, in which trajectory separation in the lifetime of an eddy, is 
not a perturbation any more. 

From Eqns. (4.3), (4.4) and (4.7), we can establish a perturbative calculation scheme for 
and v^. Notice that, within perturbation theory, is a one- valued function of x and t, and 
v(x, t) defines automatically a velocity field for the solid particles. The separation between Si and 
all the other scales of the problem, has the consequence that, in the present case, the Weinstock 
approximation is exact p4[ |. What happens is that trajectory separation is produced mainly by 
eddies of size r > S, for which trajectory separation is a perturbation. This has the consequence, 
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in particular, that the Weinstock approximation applies also at scales r < Si for which trajectory 
separation is not a perturbation at all. For dominant gravity, i.e. when uq > us, trajectory 
separation is produced mainly by the gravitational drift ug and the Weinstock approximation is 
automatically satisfied. 

We can calculate at this point the time correlation for the solid particle velocity and adopt the 
approach followed in [^5[ we can thus write, using Eqn. (4.7): 

«(0, 0)<(0, t)) = J ^^{u\^iO)u\^it) exp[ip . (z-m 0) - z(t|0, 0))]) 

^ ■ exp(ip ■ UGt) ^ ^ .. _ (4.10) 



(27r)2(27r)2 ^ " ^ Wki(0)<5Jpi(0 



j=pJt 



where 



Z[3] = (exp (il Af-y dtu^(i) . J,(t)^ 

= A/'cxp(-i I drdr'l ^J,(r) • U,-(r - r') • J-s(r')) (4.11) 
is the generating functional for the field and 

{0, T>t 
-exp(-^), 0<r<t (4.12) 

(1 - exp(-t/T5)) exp(r/T5), t < 

Substituting back into Eqn. (4.10), we obtain, after introducing dimensionless variables t — t/rs, 

UG = koTsUG and 7 = Tsjko = ts/{2tl): 

«(0,0)<(0,t)) = ^ J dx x~^MuG{x~l)^+HuG{x~l)^t)] 



X exp 



- 1 7^(2; -1) r dy / r e-^*2'= 



9 2 , V. MV- - . ^ ' (4.13) 

Jl 2/3(1+72/3)^ 1 — 72/3 

We see from this equation that decorrelation of the fluid velocity sampled by a solid particle receives 
three contributions: one from the gravitational drift ug, one from the eddy decay ^x^t and the 
integral term in the exponential, which comes from inertia produced trajectory separation. This 
last term is peculiar, in that it saturates to a constant for long t instead of continuing to increase 
indefinitely. This term is the argument in the exponential expression for Z[3] [see Eqn. (4.11)], 
which is essentially: 

pp : ([zP(t|0,0) -z(t|0,0)][z^(i|0,0) -z(t|0,0)]) (4.14) 

with the drift ug subtracted out, and with p the wavevector entering the integral of Eqn. (4.10). 
But, from Eqns. (4.7-8), we saw that this expression saturates at t ^ 00. In consequence of this, 
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for long enough times, the large x behavior of the integrand in Eqn. (4.13) will be dominated by 
the value at saturation of the inertia produced term. 



1. Velocity self-difFusion 

Inertia causes two ranges of time separations in the correlation (u^(x, 0)u^(x, t)); one at short 
times dominated by sweep from the velocity difference Uq + v — u and one at long time associated 
with eddy decay, where Eqn. (3.10) holds The transition between the two ranges occurs at 

^^max(.|,4) 

From Eqns. (4.5-6), for dominant inertia, i.e. us > uq, this cross-over time is much shorter than 
tl, while, for dominant gravity, i.e. for uq 3> us it is possible that sweep dominates for all inertial 
timescales; for this to occur, it is necessary that the crossing time of a large eddy by the particle 

3 

be less than r^, i.e. kQUQT^ > 1. For dominant inertia the cross-over time Ug/ {p^C^^it) ^ p^^Ts 
is just the lifetime of an eddy of size Si [see Eqn. (4.9)]. 

For dominant gravity, the exponential term in Eqn. (4.13) can be neglected. For t 4C 
min(TG,TL) with tq = ^(^)Vl - "3/2^ . we find: 

(K(0,t) - <(0,0)]2) = ^ y dx x-l [1 - MuGt{x - 1)^) - Hnatix - 1))^)] 

~ -a5CKo,ei(wGi)* (4-16) 

where ~ 2.149 [see Eqn. (3.25)]. The time tq, for ug < ul, is the lifetime of vortices whose 
lifetime equals the crossing time by a falling particle; for p = 0(1), tq coincides with the eddy 
turn-over time of vortices with characteristic velocity uq- 

For dominant inertia uq < us and short enough times t p~'^Ts, only the last piece in Eqn. 
(4.13) will contribute and will be quadratic in i; if ^ t^: 

e-7*!/* - e"* 1 1 ^ 



ly-'t' (4.17) 



1 — 7?/ 3 2 

Substituting into Eqn. (4.13), we are left with the following expression: 

73p(x-l) p dy 



([<(0,0-<(0,0)]2) = ^ / dxx-|{l-exp 



V J I y{l + ^y^)- 
Ax a:"3|l-exp [ ^ 42 ]} (^-l^) 
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where use has been made again, in passing, from the first to the second hne, of the condition 
t^p-^Ts. Using /o°°dx x^s (l - exp(-Aa;)) = 3r(2/3)A3 , we obtain therefore: 

(K(0,t)-<(0,0)]2) c ^T{2/3){3logiTL/Ts)y'CKoilhust)i (4.19) 

As predicted in [p2[ , at short times, the time structure function for u'' has a sub-diffusive behavior 
with exponent | both for dominant uq and dominant us- What happens is that at such short time 
scales, the particle crosses at constant speed (remember also, in the inertia dominated case, that 
Sc ^ Si) vortices whose velocity field is, in the limit, basically frozen; hence a Taylor hypothesis 
applies, and time correlations coincide with their spatial counterparts. 

2. Velocity correlation times 

Starting from Eqn. (4.13), we can calculate the correlation time Tp for the fluid velocity sampled 
by a solid particle: 

/"OO 

rp = {[ulfy^ di«(0,0)<(0,0) (4.20) 







To lowest order, any discrepancy between the PDF (probability distribution functions) for and 
can be neglected and we have ([u^]^) = {[u\]'^) — j'^t- We begin by analyzing the case of 
dominant inertia: uq = 0. Taylor expanding in p^^ the integrand in Eqn. (4.13) and substituting 
into Eqn. (4.20), leads to terms which diverge when integrated in x. This indicates that the time 
independent part of the inertia term in Eqn. (4.13) dominates the integral. We thus Taylor expand 
in p~^, only the time dependent piece of the integrand in Eqn. (4.13), i.e.: 



exp 

to obtain: 



7te 3 + ^ \ ^ ^ \ e-* + — — (4.21) 

/ 1 7/3 14- ^7/ 3)^ \ — ^7/3 /J 



V h 2/3(1-1-72/3) V i-7y 



di«(x, OX(x, t)) = -f dx a;-3 exp ( - ^^^^ ^ 



y3(H-7y3) 



'^y—TT-, IT^TT: _ 1,, 1 r- (4-22) 



-7X3 2p^ 7l y 3 (1 — 72^ 3 )(1 -|- 7x3 )(X3 -f ?/3 ) - 

and we see that the integral in x of the 0{p^^) on second line of Eqn. (4.22) is dominated in fact 
by a saddle point at a; = (k/ko)^ ^ {p/lY^ i-e- at fc ^ "^T^- Combining this result, with the fact 
that the integrands are peaked &i y ^ 1, Eqn. (4.22) will take the form: 

7^2; Ay 



I dtK(0,0)<(0,t)) = ^ / dxx-texp(-14 / . 
6 7i V 2p^ 7i „3 1 + ^„3 / 
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We thus obtain, for the deviation Tp — t^: 



= 1 + (4.24) 

TL 



where 

B{l) = (^) 'r(l/3) [i - I + f + f fog ^] ^ (4.25) 

It is to be noticed that the factor B{^) is always positive, i.e. the correlation time for the fluid 
velocity seen by the solid particle is longer than r^. Following the argument in this would 
be expected in the case of a velocity field with statistics defined in an Eulerian frame, and is 
exactly the result obtained in |38j. In the case of a Lagrangian statistics, it is not clear whether 
the deviation between solid and fluid particle trajectories, should have lead to a faster, rather than 
slower decorrelation rate. 

In the case of dominant gravity, as expected jl^, ^ , there is always a decrease of the correlation 
time. In place of Eqn. (4.22), we have: 

'OO 

di«(0,OX(0,t)) 



y2 roo i<oo 4 1 11 

= ^J dt y dx x~^3o{uG{x-l)^t) + ,h{uG{x-l)n)]exp{-^txi) (4.26) 

which, using dx3,y{(3x)e~°'^ — I3^^{a^ + /3^)^^[(a^ + /3^)^ - a]" ||3l|, leads to the expression 
for the correlation time: 

r^(-^)^[((-^)Vl)h] (4.27) 

tl o Jo x3 Vui3(a; - 1)2 / ^^^ucyx-lp^ ' 

We can obtain limiting expressions for this ratio, when the crossing time (fcouc)^^ is much longer 
or much shorter than the integral time t^: 

3. Eulerian correlations 

The limit T5 — > oo, corresponding to the case of a particle with infinite inertia, leads, from 
Eqn. (4.2), to a particle velocity, which, in the absence of gravity, is identically zero. Hence 
u''(x, <) = u(x, <) and the time statistics for the fiuid velocity seen by the particle coincides with 
the Eulerian turbulent statistics. In this regime, the dimensionless units introduced for Eqn. 
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(4.13) are not appropriate any more. Redefining t = "fkot, Eqn. (4.13) takes the form, after writing 
exp(-i/Ts) ^ 1 -i/rs: 



(ui(0,0)Mi(0,t)) 



4 r 

6 Jl 



dx X ^ exp 



- 



1 {x-l)f3i 



/oo 
dyy-^eM-iy^))\ (4-29) 



2p2 V 2 

We start by calculating the Eulerian correlation time 

TE = j di(u(x, t) ■ u(x, 0)) (4.30) 

Contrary to Eqn. (4.13), it is the linear in t, 0{p~^) term in Eqn. (4.29), which, at fixed long 
enough t, dominates for a; — »■ oo. The same reasons leading to expand Eqn. (4.21), suggest that 
we must now expand: 



exp 



(x-l) 



2p2 



(4.31) 



Instead of Eqn. (4.22), we find: 



Jo 



dt{ui{0,0)ui{0,t)) = 
(x-l) 



7*5 

6 Jl 



dx X ^ 



[x-iy^ 



/oo 
dy y~'^\x^ 



3 +y3 + 



2p2 
3(a;-l)^-l 



1 3(:r-l)^ 



4p2 



(4.32) 



2p2 7i ^ ^ ' ^ ' 4p2 

All the terms involving factors p^^ lead, after integration, to an 0{p^'^) result, except one which 
leads to a 0{p~^ logp) term; the integral in Eqn. (4.32) will read, to leading order in p: 



i: 



dx 



_ 5 
X 3 



l + /0"^a; = 



+ 0{p-^) ~ ^ - 



31ogp 



We obtain then the result for the Eulerian correlation time: 



^ = 1 



21ogp 



(4.33) 



(4.34) 



TL P 

which is shorter than tl, as expected from the fact that the velocity field statistics is defined along 
fluid trajectories, and, sampling at fixed space position should lead to an increase in the rate of 
decorrelation. Comparing Eqns. (4.28) and (4.36), we see therefore that there is a transition from 
a correlation time longer than tl for light particles, to a shorter one for heavy particles. The origin 
of this lies in the opposite orderings rs < t and ts ^ t, on which the Taylor expansions of Eqns. 
(4.21) and (4.31) are based. [More precisely, for rg > ptl, we have koSi > 1 and the saddle point 
in Eqn. (4.22) disappears]. 

As a last exercise, it is possible to calculate the sweep produced decay in an Eulerian two-point 
two-time structure function in the form: 



Srr{r,t) = {[Ur{r,t) - Ur{0,t)][Ur{r,0) - Ur{Q,0)]) 



(4.35) 
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From the discussion leading from Eqn. (3.14) to Eqn. (3.17), one finds that the structure function 
in Eqn. (4.35) is obtained by inserting a factor 2[1 — Joi'f'x) ~ J2(?'a;)] in the integrand of Eqn. 
(4.29). If one considers shorter time and space scales kor <C 1, t ^ tl, the leading cause of 
correlation decay is sweep, and the tx^ in the integrand of Eqn, (4.29) can be disregarded. Again 
because of shortness of i/r^, one can Taylor expand exp(— ty^) in the same equation and the final 
result is: 

Srrir,t) ^2CKoieir^ dx x"! [1 - Jo(x) - J2(a;)] exp - ^ „ (4.36) 



The term in the exponent is ©(i/T^-i)^, with Tj,-i the sweep time at scale r. Hence, if < 3> T^-i, 
it is possible to Taylor expand the Bcssel functions and the result is 

Srrir,t) - Srr{r,0) J dx x^ cxp ( - {t/Tr-i)^X^j - S'„(r, 0)(^j ' (4.37) 

i.e. a power law decay of the structure function for times longer than the sweep time at that space 
separation. 

V. Solid tracers: concentration fluctuations 

Because of inertia, the particle velocity field v(x, i), contrary to u(x, <), does not preserve 
volume. Physical intuition suggests that particles which are denser than the fluid, will tend to 
concentrate near the instantaneous hyperbolic points of the flow, and to escape from the elliptic 
ones |3^. For this reason, a distribution 9{x,t) of solid particles, in the absence of external 
sources, will be characterized by finite amplitude fluctuations superimposed to a uniform mean 
concentration field 6. These fluctuations are expected to have a correlation time of the order of ts 
and a correlation length determined in consequence. We are going to neglect any effect of gravity 
and set from the start ug = 0. We will also limit our analysis to the case in which ts is in the 
turbulent inertial range, i.e. we consider ts tl (more precisely, ts < p~'^tl). In this way, all 
non-universal effects associated with the large scales of the flow are eliminated from the problem. 

The length Si is crucial to the two-particle statistics, in that it gives the scale below which 
solid particles move ballistically relative to one another. In fact, Sc fixes the cross-over scale to 
ballistic behavior, only for the relative motion of solid and fluid particles; the resulting picture 
is given by pairs of particles, separated by Si, moving ballistically over scale Sc- It is easy to 
see this: if A,.t; is the typical relative velocity between two solid particles at separation r and 
ArU ^ ^Koii^^)^ ^^'^ corresponding value for the fluid velocity, one will have for r <^ S, from 
Eqn. (4.2): A^u ^ (rs7r-i)~2 A^u; exploiting the fact that the characteristic time of variation for 
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V is TS, the condition ts^j-v ~ r, gives then r ^ Si. 

The concentration correlation ©(r) — {9(r,t)9{0,t)) is proportional to the equilibrium PDF 
P(r) for the separation of a pair of solid particles advected by u(x, f). The separation r{t) obeys 
an equation in the form r{t) = v^(x + r, t) — v^(x,t) [we use from now on the shorthand r(t) = 
5z^{t\x + r, 0)], and, for r ^ Si, the separation process takes a diffusive nature: 

^([r„(i) - r^mM) - MO)]) = 2DMr) (5-1) 

A finite level of concentration fluctuations, in the absence of external sources, is associated with a 
finite divergence of the diffusivity tensor: dcDa^ ^ 0. If this component of the diffusivity tensor is 
small, it is possible to proceed perturbatively: = D'~°}+D'~^}, P = P""+P''', with 0^0^°''^ = 0, 
P(°> uniform and P<i'(r) a ([6»(r,i) - 0{O,t]'^); the equation for the fluctuation amplitude P<-^^{r) 
would read therefore: 

^i°^5„5^P'^> = -P'^'dadpD'^j, (5.2) 

The procedure to determine is similar to the one leading to Eqn. (3.17). Prom Eqn. (4.4) 
and the relation r{t) = v^(x + r,t) — v^{x,t), we obtain: 

D„,(r)= lim i ["dU fdt. r ^ r ^exp(- *^+^^-^^-^^ )5-,(r,n,r.) (5.3) 

T^oo T Jo Jo J_^ TS J-oo TS V TS ' 

with the time correlation of velocity differences along solid particle trajectories: 
5^^(r, tut-i) = (K(r, h) - <(0, h)] [u}ir, t^) - u}iO, ts)]) 

= 2[{ul{r,h)u}{r,t2)) - {ul{r,t,)u}{0,h))] (5.4) 

We notice that, if we approximated S'^^(r, ti,t2) = S^^{r, ti — 12), since daS^f^{r, ti — 12) = 0, we 
would obtain from Eqn. (5.3) a divergenceless Dai3{r). We have to take into account therefore the 
effect of trajectory separation described in Eqn. (4.7). Proceeding as in the case of the 1-particle 
statistics, we arrive at the following modification of Eqn. (4.10): 

{u^{O,ti)u0{r,t2)) = - j^-pr-r^eM^p-r) 5.5 

where 

Js,rHt, (r) = Jt, (r) - e'^ Vt, (r) (5.6) 

and Z[3] and Jt are given in Eqns. (4.10-11). Carrying out the wavevector and time integrations 
in the definition of Z[3] and using Eqns. (5.6) and (3.1) leads, after some algebra, to the following 
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expression for the velocity correlation: 



Jo 



L 



2tt 

/o 



cos^ (j) + (SajS — — y-) sin^ (j) exp(iA;r cos (/>) 



X expj - CKoi^^rlk^ j 



[F{s,ti,t2)+ G{s, tiM) ( Jo (sr) + J2 (sr) cos 2(^)]sds j 



(fc2 + s2)f(l + 7,rs) 



where is the angle between k and r, 



(5.7) 



F(s, tiM) = f{s, h) + f{s, t2), G{s, tiM) = f{s, h - t2) - F{s, tiM) 



(5.8) 



and 



f{s,t) 



-\t\/rs 



-\t\/rs 



IS 



(5.9) 



The effect of trajectory separation is contained in the last line of Eqn. (5.7). We see that the 
contribution, which leads to finite divergence of the correlation (m^(0, ii)u^(r, 12)), is the (f) depen- 
dence of this factor. The remaining (j) dependence, contained in the second line of this equation, is 
simply the factor kj^kj^ exp(ik • r) arising in the Fourier transform of Eqn. (3.1), and would give 
by itself zero divergence. 

The argument of the exponential in the last line of Eqn. (5.7), for fixed ts/tl, is 0{p~'^), so 
that we may try a Taylor expansion. However, as it happened with Eqns. (4.21) and (4.31), the 
resulting integrals in k diverge. Wo therefore keep in the exponential the leading contribution in 
k, which is the time independent piece of its argument, and expand the remnant, which, to leading 
order in p, gives the following expression: 



1 — CKoi^^Tgk^ cos 2^ exp | 



1 — Cxoi^^Tgk'^ cos2(/)exp | — 



(/eg + s2) 3(1 + 7, rs)J Jo (/eg + s2) 3(1+ 7, rs) 
SCKoieiji^X f°° G{s,ti,t2)h{sr)sds 



2k^ 'Jo (fcg + s2) 3(1+ 7^ Ts) 



(5.10) 



plus terms which would lead to a divergence free contribution to (^^(0, ti)u^(r, f2)) and would 
disappear from Eqn. (5.2). Substituting into Eqn. (5.7) and then back into Eqns. (5.4) and (5.3), 
we find, after carrying out the time integrals and the integral in </>: 



r)(i) 



f 

Jo 



X 3 dx exp • 



2r2 



dyy =*J2(y) 



1 



2(l + (x/2/)f) 



1 \ V Tr 



(5.11) 
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and it is possible to see that -D^"^ is given by the same expression vahd for a fluid parcel, i.e. by 
Eqn. (3.17): 

z,»„.lf!i^4[!iI£ + I(*,,-I^)] (5^12) 



The physical content of the expansion leading to Eqns. (5.11-12) can be clarified, noticing that, in 
a way perfectly analogous to Eqns. (4.10-11), the generating functional Z[3] entering Eqn. (5.5) 
can be written as: 

(exp (ik- J dt(u^(0,t)4(t)-u^(r,t)4(t)))) 

^exp{-^ : I dtdt' [U(0,t-t')(4W4(i') + 4W4(i')) + 2U(r,t-i')^tiW4(i')]} 

The argument in the exponential is in the form A:^C/(r, 0)r| ^ k'^{Sf + CKoii^r)^ Tg), where the 
term involving Si gives the one-particle contribution to trajectory separation, while the remnant is 
the two-particle correction coming from r > 0. Substituting into the definition of 0^/3 gives then, 
using r Si: 

J kdkUkJk'e^'''^' [l - exp ( - k^CKoi{er)iT^s - ik ■ r) 
Mfc[/fc7-'(l-e-"^"-)+ / k^dkUk%'CKoi{er)iT^e-^'''" 



which, using Eqn. (4.9), is ^ — [ra -|- p S^ {-^^)^]- Thus, the origin of the Z?*^' in the two- 
particle contribution to trajectory separation is confirmed, together with its being generated at the 
/cQ-dependent scale Si. This is of course confirmed by direct analysis of Eqn. (5.11); the integral 
is dominated hy k = x/r ^ Sj^^ and we obtain, for r > Sf. 

D^'>^ {rSi)i ^ ^tlJ ^ ' 

Thus, for Ts < p^'^TL, D^'> < D"" and Eqn. (5.3) applies. 

Using the relation: J^'^ dxx-h'"""^ 34{x) = ^g^^^^^M{7 -1/ (Aa)), with M{a,b,x) the 
confluent hypergeometric function pO[ , we obtain in general, from Eqn. (5.11): 

D'^l,^AparCl^^r^sf['^D{r/Si) + 6^pDir/Si)), (5.14) 

where we can write: 
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with c ~ 1 for r ^ Si and c ~ ^ for r <C S*; ; £) will be shown not to contribute to the concentration 
correlations. 

We can now calculate the probability P<^>(r). Substituting Eqns. (5.13-15) into Eqn. (5.2), 
after a few manipulations, leads to: 

drf'WfP''' = -p'fP'"' ' {dl{b{v) + b{f)) + -df{2D{f) + D(f))) (5.16) 

where f = r/5/. Hence, for Si <^r Si: 

P<i)(r) = p2p(„, I dyy-i [b{oo) - ydy{b{y) + b{y)) - b{y) 



r/Si 



^yp^o^bioo)[^y (5.17) 



Using Eqn. (5.15) and the limiting form for the confluent hypergeometric function M(a,6, — z) = 
P^^z-'^(1 + O(z-i)) @, we get the final resuh: 

Qir)=9^[l+f3p'[^y) (5.18) 

where 

P=4^£^^2.U (5.19) 
2la7r(8/3) 

In conclusion, we have a range of separations Si ^ r Si, in which the fluctuation correlation 
grows with a power — |, to reach amplitude ~ at r ~ Si. 

The picture which arises is one of concentration fluctuations produced at scale Si, by com- 
pressibility of the solid particle flow, and then transported to small scales and amplified by the 
incompressible part of the flow. The process is different from that of a passive scalar forced at 
large scale, due to the derivatives in the source term [and in fact the scaling exponent is different; 
compare with Eqn. (3.36)]. This source term is basically V^-D*^'(r), with Z?<^'(r) saturating at a 

4 

constant for r ^ Si and going to zero in the opposite limit. From here, the r~3 scaling of Eqn. 
(5.18) arises by dimensional analysis. 

What happens when r ^ Si? At such short distances, the separation process is ballistic and we 
cannot use a diffusive approximation anymore. In p^ , it is suggested that the correlation build-up 
should stop only because of discreteness effects or because of the Brownian motion of the solid 
particle. Actually, extrapolating the results of the present paper to the real turbulence regime 
p — 0(1), there is good reason to think that, for ts > t^j, this build-up could stop much earlier, 
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and precisely at r ~ 5^, which, for p = 0(1), coincides with the size of vortices with eddy turnover 
time Ts- 

At separations below Si, Eqn. (5.18) ceases to be valid, and full analysis of the distribution 
P(r, Aj-v) is needed. A singularity of P{r) at r = would require focusing of A^v along r for 
r <C Si] the mechanism is sketched in Fig. 1. This means that P(r, A^v) itself should develop, 



Figure 1: Sketch of the behavior of particle relative velocities inside a domain of size R < Si. 
Particle 2 moves with respect to particle 1 at constant velocity. In order for the particle density 
to diverge as r — > 0, it is necessary that the distribution of velocities be peaked too at ^ = 0. 

as r ^ 0, a singularity at 6, where 6 is the angle between A^-v and r. The necessary trajectory 
focusing can be produced only by the compressible part of v. However, for r < Si, the production 
term for the compressible part of v can be estimated directly from the second term in Eqn. (5.10), 
to be 0(p~^) relative to the rest, and is able to act only for a time ts in the ballistic region. Hence, 
P(r, v) must be singular before this region is reached. On the other hand, for r > Si, where the 
diffusive approximation works, the asymmetry of P(r,v), associated with compressibility of the 
flow, can be estimated from |^°g|^|(o) ~ D^^^ /D^°^ so that, if rg < p~^tl, singularities should not 
be expected in P(r, v), for ^ = and r > Si either. The conclusion is that a plateau for 6(r) 
should be present at r < 5,. 



One of the consequences of the compressibility of the velocity field v(x, t) is that the ergodic 
property is not satisfied anymore: velocity moments calculated along solid particle trajectories 
differ from those obtained from spatial averages. As mentioned before, physical intuition suggests 
that solid particles should privilege in their motion certain regions of the fiuid with respect to the 




2 
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others (namely, hyperbolic with respect to elliptic regions). It is difficult, however, to translate 
this into a statement on the form of the PDF for the velocity u'^. 

We have at our disposal the equations satisfied by the velocity field u''. It is possible therefore 
to calculate its moments and to reconstruct its PDF. We consider the case of zero gravity Uq = 
and ts/t]^ small. As in the analysis of the concentration fluctuations, all non- universal effects 
associated with the large scales of the flow are thus eliminated from the problem. From definition 
of u'' and Eqns. (4.1-2), we obtain the following set of equations, valid to lowest order in p^^: 

i {dt + u(x, t) ■ V)u-(x, t) + J d2/2^(x ~ y)u-(y, t) ^ J d^yhi^ - y)^(y, t) 
\u(x,t) ^ u-(x,i) - v-(x,i) - /!^dTexp(-^)u-(x,t) 

which differs from the analogous equation for because of the non volume-preserving advection 

term u ■ Vu'' . From here we can carry on standard field theoretical perturbation theory, either by 

the Martin-Siggia-Rose formalism ||2^, or working directly with Eqn. (6.1). The building blocks 

of the diagrammatic expansion are shown in Fig. 2, and are the propagator Gka/s: 

/■* k kg 

Gka/jUk/sW / dr exp(-7fc(t-r))-|2-Uk/3('r) (6.2) 



the correlator U^^pi^) 

(fc2+fc2)f 

and the vertex Fka/j^: 



^kaffW ^ ^r^ .,, , ,,.4 exp(-7fc|t|) (6.3) 



Ti^apjUpi^ul^it) ^iXspSa-ySik + p + s) dr exp( )ul^{T)drUpp{T) (6.4) 

J — oo '^S 

where the coefficient A 1 is introduced, as in Eqn. (2.15), only for the purpose of book-keeping. 
To lowest order in A, the correlations for the fields u^(x, t) and u''(x, i) are trivially equal. To 
higher orders, differences arise, which would not lead, if V ■ u = 0, to differences between the one- 
point PDF's for and (see also ^^). In our case, this is not so, and the difference between 
the moments of the two PDF's can be calculated in perturbation theory; to 0(A"): 

(K)")- - / (6-5) 

where (uki (^)---i*k,„ (i))'"' is the sum of the Feynman diagrams with m outgoing velocity lines and 
n vertices. Because of symmetry under space reflection, the lowest order contributions are ©(A^); 
the corresponding diagrams are shown in Fig. 3 and lead to corrections to the velocity second 
and fourth moments. In order to check for the presence of divergences in loop diagrams, we carry 
on power counting on Eqn. (6.1). Rescaling space and time as in Eqn. (2.16), we find [A] = 0, 
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a 



b 
c 




Figure 2: Fcynman diagrams for the propagator Gka/3 (a), for the correlator U^^^ (b) and the 
vertex rka/37 i'^)- 




Figure 3: Feynman diagrams providing the fowest order correction to ((w^)^) (a) and {{u'')'^) (b). 
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implying the possibility of logarithmic divergences. Now, a perturbation expansion in A of Eqn. 
(6.1) ceases to be sensible at scales below the length Si defined in Eqn. (4.9). Prom Eqn. (4.13), 
the effective decay rate for the field appears to be: 

P ilk kSi^l 
= [usk kSi » 1 ^'-'^ 

and this expression should be substituted for 7 in Eqns. (6.1-3). For kSi » 1, 7^^ is just the 
inverse of the crossing time of an eddy of size k~^. In this large k range, the appropriate scaling 
for the frequency should be, instead of the one provided by Eqn. (2.17), which lead to [A] = 0, the 
following one: 

[t] = l, [X] = -[u] = -^ (6.7) 

The change of scaling in 7'' is therefore sufficient to regularize the divergent diagrams, providing 
an effective ultraviolet cutoff at k = S^^. 

We calculate explicitly the loop diagram in Fig. 3, and the corresponding correction to {{u^)^): 

= / 7?^7^ /° dii /° dt, r dri t 



(27r)2 (27r)2 

X exp - 7^(ii + t2) - ^J'^ —'^S„0dr^dr^SjSs[U^^siTi - T2)Ul„0{ti - tl) 



where k = — p — s. The time integrations can be carried out at once and, after some algebra, we 
reach the following result: 



P\2\(2) 



, ( d^p d^s 

7 



d2p d2_S £^;^j7g(P-L-s)' 

0' (7j^ + % + 7?)(7^ + 7? + Ts^){% + Ts^hlv" 



s 7? + T-s -I 
where Uk = C^kaa(O). As predicted in the discussion leading to Eqns. (6.6-7), substituting 7"° — > 7 
would lead to a logarithmically divergent integral. Comparing with Eqn. (4.9), we see that this 
integral receives contribution from wavcvectors in the range [S~'^, S~^], i.e. from those eddies fast 
enough for the particles to be unable to respond to their velocity field, but still sufficiently slow 
for trajectory separation to be considered a perturbation. To find the leading behavior in S^^, the 
integral can be rewritten, after the change of variables y = (7 , 1 r^)"^, z = p/q, in the form: 



167r3p2 ^ J^_^ y + si + fci)(fci +si + y){pi + y)kipi 
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/ 2 2 -2 . yCOSm, 2-2 si On , n 

-(^3+s3 +fc3+y) + ^^^ — "Lis^+k^-y) +0{p-^) (6.9) 

S3 + y J 

where k = (ps)^^k, p = (ps)~'p, s = (f>s)~^s, coscj) = p • s. We obtain then the final result: 

{{un'r = ^4 (6.10) 

where: 

3 r sin^ 6 



Stt^ io io ^ (2; + 2;-! + 2cos(?i)H {z + z-'^ + 2cos(l))i + zi 

cos(h 1 
— r ^-0.32 (6.11) 

Z3 + z 3 + (^z + z~^ + 2cos(j))3 ^ 

is evaluated by numerical integration. The correction to the velocity amplitude is negative. In the 
presence of inertia, solid tracers prefer therefore to lie in regions of the flow where the turbulent 
velocity is smaller. 

Extrapolating Eqn. (6.10) to p = 0(1) suggests that {{u^)"^} — ^ w|. We can have some 
idea of what we should expect for dominant gravity us < uq from dimensional analysis of Eqn. 
(6.8). In this case 7^ uck, the inverse sweep time due to the particle fall, and we would find 
((u'°)^)'^' ~ with ~ ugts giving the transition to the small scales for which the sweep 
time is shorter than rg, and to which the particles are unable to respond. Prom here we find 
((m'°)^) — ~ {us/uG)^u'g and we see that gravity reduces the amount of non-ergodicity of the 
solid particle flow. 

VII. Conclusions 

Consideration of a finite correlation time in the transport by a random velocity field has allowed 

analysis of a series of issues. We summarize the main results: 

i The self-diffusion of a fluid parcel obeys linear scaling in the inertial range (as it should) 
with a universal constant Cq = C^^i^^logp/ p [see Eqns. (3.7-3.10)], which is sen- 
sitive both to the ratio of the eddy turn-over and life time, and to the rate of eddy 
velocity decorrelation at times much shorter than the eddy lifetime. A quadratic max- 
imum (at least), at time separation equal to zero, is necessary for Co to remain flnite. 
(An exponential time correlation, for instance, would not satisfy this condition). This 
sensitivity on the short time behavior of time correlation was not observed in any of 
the other transport processes considered in the present paper. 

ii The relative diffusion of a pair of fluid parcels, exhibits (again as it should) Richardson 
and normal diffusion behavior, respectively, for coordinates and velocities. The PDF 
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for relative separation is a stretched exponential with exponent | [see Eqn. (3.21)] and 
it is possible to express the universal constants c and c entering respectively coordinate 
and velocity dispersion, in terms of the parameter p. Precisely: c ~ °' "^^^ and 

3 

c ~ S-OSTCj^^i jgg^ Eqns. (3.23) and (3.26)]. For the Batchelor constant, we obtain 
instead [see Eqn. (3.39)]: B ~ 11.32p. 

iii The correlation time rp for the fluid velocity sampled by a solid particle has a behavior 
consistent with previous analysis neglecting the structure of the turbulent inertial range 

|3^ . Values of tp/tl, above unity are found for dominant inertia and ts ^ tl, with 
tp/tl — 1 = [see Eqn. (4.25)]. On the contrary, in the case of dominant 

gravity, rp /r^ < 1 irrespective of the value of the ratio ut / uq between the turbulent 
and the fall velocity; specifically [see Eqn. (4.28)], we find tp/tl — 1 = 0{{uGkoTL)'^) 
for ug <IC ul, and tp/tl — ©((ug/coTl)^^) in the opposite case. For short times, the 
expected sub-linear behavior for the fluid velocity along a solid particle trajectory is 
found: {\u^{x, t) — u^(x, 0)|^} ^ {luAt)i ^ with A = G^S depending on whether gravity 
or inertia dominates [see Eqns. (4.16) and (4.19)]. 

iv The Eulerian correlation time te (and by continuity, therefore, also Tp, in the regime 
Ts 3> T]^) is shorter than its Lagrangian counterpart, with teItl = 1 ^ 2p^^ log/O [see 
Eqn. (4.34)]. Sweep produces a power law decay of correlations between velocity incre- 
ments in the form Srr{r,t) — {[ur{r,t) — Ur{0,t)][ur{r,0) — 1*^(0,0)]). More precisely, 
for time separations longer than the sweep time T^-i: Srr{''',t) ^ Srr{''',0){Tj.-i /t)^ 
[see Eqn. (4.37)]. 

V In the absence of gravity, and for p^r,, <^ rg <C p^^tl, the spectrum of concentration 
correlation induced by turbulence in a solid particle suspension, is universal and has 
power law behavior for separations above the size Si of an eddy which is crossed by a 
typical solid particle in a time equal to its lifetime. More precisely: 6^"^ {6{r)6{{))) — 1 ~ 
/92(5,/r)i [see Eqn. (5.18)]. 

vi The solid particle flow is non-ergodic, with a difference between the fluid velocity sam- 
pled along a solid trajectory and the corresponding Eulerian average: {{vFY) — ~ 
_oj32k3gp^2 jggg Eqns. (6.10) and (6.11)]. Dimensional reasoning for p = 0(1) sug- 
gests that gravity should reduce this effect from {{u^)'^) ~ ^ u'g to {{u^)'^) — v3p ^ 
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(m5/"g)'us- 

Analysis of some of these problems, actually did not exploit the finite correlation time of the 
velocity field produced through Eqn. (2.7). In particular, the process of fiuid parcel relative 
dispersion was considered to the same order in p as in the Kraichnan model and, to this level, no 
information on the Lagrangian statistics was necessary. Finiteness of the correlation time had the 
only purpose to allow a meaningful definition of quantities such as CkoI and e. 

In the case of the self-diffusion properties of fluid and solid particles, a finite correlation time and 
inclusion of the Lagrangian nature of time correlation was necessary from the start. Nonetheless, 
the only point in which analysis of the random velocity field could not be avoided, was to determine 
the dimensionless constant Cq jlj, |2^; the diffusion exponents in the various cases were already 
available by dimensional reasoning. 

Evaluation of the correlation time rp and analysis of concentration fiuctuations and non- 
ergodicity of particle trajectories (points iv-vi), instead, rested heavily on the fact that the cor- 
relation time was finite and on knowledge of the actual form of the random velocity field time 
correlation. The analysis confirmed the role of eddies with lifetime rs, already pointed out in [ p2| . 

Some comments are due on these last issues. As regards correlation times, they depend in gen- 
eral on non-universal aspects of the velocity statistics, and, in the present case, on the assumption 
that also the large scale statistics is defined along Lagrangian trajectories. In consequence of this, 
the Eulerian time of the fiow resulted shorter than the Lagrangian correlation time. (Following 
prf , the Eulerian correlation feels, at the same time, the decorrelation from relative motion of 
the fluid, and the effect of eddy decay). For ts ^ t^, the standard picture of inertia and gravity 
leading, respectively, to increase and decrease of the correlation time, however, was confirmed. 

As regards concentration fluctuations, previous treatments of this problem, either were limited 
to the case of particles with Stokes time shorter than the Kolmogorov time of the fiow |24 , or 



neglected turbulent small scale structures altogether ||2^ . This was due to the difficulty in analyzing 
trajectory crossing effects on inertial range scales, associated with the need for a proper treatment 
of the Lagrangian time statistics. The fully kinetic treatment adopted here, in which the relative 
motion of individual solid particles is fully taken into account, in contrast with the fiuid equation 
approach used in ||2^ , together with the large p limit, is what allows treatment of the problem. 

It should be mentioned that solid particle concentration fiuctuations may be important in the 
process of rain formation. It is known that the settling rate of a suspension is enhanced in the 
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presence of clumping of the heavy particles , and turbulence induced concentration fluctuations 
appear to be one of the important actors in the process |43j. Inclusion of the effect of gravity, on 
the same lines of the analysis carried on in section IV would therefore be necessary. 

As regards non-ergodicity of the solid particle flow, it should be mentioned that this is a problem 
one has to deal with, before trying to extend standard Lagrangian transport models (in particular 
the well mixedness hypothesis on which they are based |]l^) to the case of solid particles. 

An important aspect that must be stressed, in the calculation of both rp and the concentration 
correlation spectrum, is the role played by the localization length Si . This length ceases to have a 
physical meaning for finite p, nonetheless, it fixes, in perturbation theory, the scale at which both 
fluctuations and the difference rp — tl are generated. Notice that, in the case of concentration 
fluctuations, this occurs in spite of the fact that the concentration correlations are peaked at the 
inertial scale Si. 

Another peculiarity of the large p expansion is the multiplicity of space scales associated with 
eddies having time or velocity scales related to ts and U5 [see Eqns. (4.6) and (4.9)]. All of them 
collapse, for p — 0(1), on the size of a vortex with turnover time equal to ts- In real high Reynolds 
number turbulence, this is the saturation length expected for concentration fluctuation build-up, 
when Ts is an inertial range quantity. 

The parameters p and p are central to the extension of the Kraichnan model to finite correlation 
times. The situation of reference in real fiows is the inverse cascade range of two-dimensional 
turbulence. An estimate of these parameters could be obtained using the leading p expressions 
provided by Eqns. (3.10), (3.23), (3.26) and (3.39), with the values of the constants Co, c, c and 
B obtained from DNS. For instance, assuming p = 00, comparison with the results presented in 
H would give p ~ 2. 

The results of the present paper have been obtained to leading order in p. To this order, no 
perturbative effects in the structure of random velocity fields are present, and the correlations 
for the Lagrangian velocity obey Eqn. (2.6). The parameters entering these correlations must 
nonetheless be considered as renormalized quantities in a renormalized statistical field theory. No 
claim on the nature of these renormalizations is made, apart that, to lowest order, marginality of 
interactions suggests that correction to scaling be only logarithmic. 

To this order in p, extension of the results to three dimensions presents no conceptual difficulties. 
In particular, the mechanism of production for concentration fluctuations, and for correlation time 
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and PDF corrections, is not expected to suffer modifications. Whetlier a random velocity field 
model like the present one, could be appropriate to describe transport by a three dimensional 
turbulence, laden with coherent structures and intermittency, is a different matter. 

The present extension to finite correlation times of the Kraichnan model is perturbative in 
nature. Imposition of time statistics along Lagrangian trajectories had as consequence a non- 
Gaussian velocity field. This resulted in a field theoretical perturbation theory, with expansion 
parameter p~^, which is somewhat different from other field theories arising from closure analysis of 
the Navier Stokes equation. It would be interesting to understand the relation with such theories, 
in particular with the quasi-Lagrangian approach described in and following papers based on 
this work (see and references therein). 

There are situations in which the higher orders in p^^ become necessary. A relevant example 
could be the derivation of a turbulent closure: in this case, extension of the theory to realistic 
values of p could not be avoided. Related to this issue, is the calculation of the anomalous scaling 
exponents for a passive scalar advected by a random velocity field with finite correlation time. The 
analysis of pair diffusion carried on in section III proceeded, at the end, as if the velocity field 
had zero correlation time. To lowest order in p~^, the same zero-mode structure of the Kraichnan 
model is therefore expected To proceed in a consistent way, one should go to higher order, at 
the same time, in the passive tracer part of the problem and in the field theory for the velocity 
field. Such issues, concerning the nature of the field theoretical perturbation expansion, will be 
analyzed in a separate publication. 
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